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The functional form of open Wilson lines 
in gauge theories coupled to matter 
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The open Wilson lines are gauge-invariant operators made with a gauge transporter along an open path satu- 
rated at the end-points with matter fields. Here it is shown that numerical experiments on 3D Z2 Higgs model 
provide useful guidance in addressing the problem of the functional form of their vacuum expectation values. It 
turns out they satisfy, as long as their size does not exceed the string breaking scale, a remarkable factorization, 
related to finite-size scaling of the Fredenhagen-Marcu order parameter. This leads to conjecture a decay area law, 
like ordinary Wilson loops, with the difference that the boundary conditions of the confining string world-sheet 
are fixed along the path of the gauge transporter whereas are free on the straight line joining the two end-points. 
A consistency condition fixes uniquely the contribution of the free boundary, which turns out to be proportional 
to the string-breaking scale. Universal shape effects produced by quantum string fluctuations are also studied. 



1. INTRODUCTION 

The most basic observables of gauge theories 
are the Wilson loops. These are operators W 7 as- 
sociated to any closed path 7 of the space-time. 
Their importance stems from the fact that they 
serve as order parameters of confinement. The 
confining phase shows up in an area law for the 
vacuum expectation value of large Wilson loops: 
if 7 is scaled up keeping its shape fixed and in- 
creasing the encircled minimal area A, then (W y ) 
vanishes exponentially with A: (W 7 ) oc e~ crA , 
where the physical quantity a can be interpreted 
as the tension of the confining string where the 
chromoclcctric flux is concentrated. When dy- 
namical matter is added to such a system this 
string breaks, because pairs of matter particles 
form and act as new end points of the confining 
string. The broken string state describes a bound 
state of a static colour source -the fixed end of the 
string- and a dynamical matter field -the free end 
of the string. 

Contrary to earlier expectations, the string 
breaks only when its length exceeds a certain 
threshold value TZb- The reason, as we shall see 
shortly, is very simple. At intermediate distances 
the behaviour of the Wilson loop is indistinguish- 
able from that of the confining phase of the pure 



gauge theory. Nonetheless, when the Wilson loop 
is very large it should decay with a perimeter law, 
hence one can consequently assume that 



(W^-Cue-'^-^+Cbe-^ (1) 

where I7I is the length of 7. The first term, cor- 
responding to the unbroken string regime, should 
be implemented with a shape-dependent factor 
due to quantum fluctuations of the confining 
string Such a universal shape effect has been 
observed also in a gauge theory coupled to mat- 
ter |2] as a further support to the conjecture that 
the confining string in this more general medium 
behaves exactly like the string in the pure gauge 
vacuum as long as it is shorter than TZf, . Eq.Q 
can be used to describe the main features of the 
string breaking mechanism as a level crossing phe- 
nomenon. In particular the string breaking scale 
turns out to be 

R^2^ , (2) 



as expected [5]. Notice that /x and A, being self- 
energy contributions of the static sources, are UV 
divergent, however these divergences cancel in 
their difference, hence TZb is a meaningful physical 
scale even in the continuum limit. 

From a computational point of view it is very 
hard to observe string breaking directly through 
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the estimate of large Wilson loops, because Eq.QJ 
implies @J?] that in the loops of finite size such 
a phenomenon is visible at distances much larger 
than where the signal is drowned in the noise. 
As a matter of fact direct observation of string 
breaking by fitting the numerical data to Eq.Q 
as been reported only in two 3d systems: the 
SU(2) gauge theory with adjoint sources j^j and 
the Z 2 Higgs model g] (see Fig[TJ. 

In gauge theories coupled to matter one can 
construct gauge invariant operators that are more 
general than the Wilson loops. They are formed 
by a gauge transporter along an open path 7, sat- 
urated at the end points by the matter fields. We 
call these operators open Wilson lines. An impor- 
tant example is the U-shaped operator IA (r, t) = 

j, where the black dots represent the matter 
fields. It can be considered as an operator con- 
necting the string state — of length r at the bot- 
tom t = with the broken string state • • at the 
top, hence it is expected to have a large overlap 
with the broken string state. This suggested en- 
larging the basis of operators to extract the static 
potential by considering the correlation matrix 



3-D I2 Higgs model 



C(r,t) = 



( 



\ 



J 



(3) 



In this way a rather abrupt crossover between 
string-like and broken string states has been 
clearly seen at the expected distance IZb- So far 
string breaking has been verified in various theo- 
ries coupled to a fundamental scalar field in 2+1 
[8] and 3+1 dimensions UJ, in QCD with two 
flavours [TU|, as well as breaking of the adjoint 
string in 2+1 ^21 an d 3+1 dimensions |13| . 

It is natural to ask what is the functional form 
of open Wilson lines which generalises Eq.JIJ. 
Experience suggests that the answer to this ques- 
tion should not depend on the gauge group nor 
on the type of dynamical matter. Thus I per- 
formed an exploratory study on a particularly 
simple system, the Z 2 Higgs model in three di- 
mensions, probing open Wilson lines of various 
shapes. It turns out they fulfil remarkable fac- 
torization properties (see Eq.s © and J3J below) 
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Figure 1. Observing string breaking with a Wil- 
son loop in 3D Z-2 Higgs model. The steepest 
line is the area term of Eq.JIJ; the other is the 
perimeter term. 



They lead to conjecture a generalisation of Eq.Q 
which can be applied to any large planar open 
path 7 joining two points placed at a distance r 



(G r ( 7 )) ~ Cie -°^-\\-y\-K b ar/2_ 
pci e ~° A i- X h\+K b <rr/2_ 



(4) 



where 7 = 7 + I is the closed path obtained 
by joining the two end points of 7 with the 
straight line £, A y is the area of the planar sur- 
face it encircles, p is (a generalisation of) the 
Frcdcnhagcn-Marcu order parameter j2] , and the 
ellipses denote sub-dominant terms at intermedi- 
ate distances which become visible only at large 
scale, where the asymptotic behaviour is domi- 
nated by the perimeter law, like in Eq.Q. The 
first term of Eq. J3J is similar to the first term of 
an older proposal [Hj with the important speci- 
fication that the coefficient of r in the exponent, 
which can be regarded as the self-energy contribu- 
tion of the constituent Higgs particle, is here di- 
rectly related to the string breaking scale. When 
7 = £ then = and we get the straight line 

operator • • = G r (£) = C r . In such a case 

Eq. (gj is replaced by the following asymptotic ex- 



3 



r 




Figure 2. Open Wilson line. The solid dots repre- 
sent the matter fields; the shaded region denotes 
the world-sheet of the underlying confining string. 



pansion 

(C r )~ Cl e-» r + ■■■ (5) 

where \x is, like in Eq.JQ), the mass of the ground 
state of the so called static-light meson. Eq.(QJ 
can be still interpreted in terms of the underly- 
ing confining string (see Fig|2J). Clearly this has 
fixed boundary conditons only along 7, while it 
is free to vibrate along I. Also in this case we 
expect that the area term should be modified by 
universal shape effects due to quantum fluctua- 
tions of the string. An explicit form of such a 
contribution in the case of the U-shaped operator 
is calculated in § 0J 

2. ALMOST CLOSED WILSON LINES 

The Z2 Higgs model which has been analysed 
in detail by a number of Monte Carlo simulations 
is defined through the action 

3 

S==-p G J2 U r-frl212 c t ) * U ^<t>x+li • (6) 

P X l_L=l 

Here U x ^ and Up are link and plaquette variables 
respectively, which describe the gauge degrees of 



freedom; the site variables <p x represent the mat- 
ter fields. Both U x p, and <f> x take values on the 
Z2 gauge group. This model is self-dual and has 
a confinemcnt/Higgs phase and a Coulomb-like 
phase ^3] . All the MC simulations are performed 
in the region < (i G < 0.755, < fa < 0.248 
which is called confinement "phase" , because the 
Wilson loops of intermediate size decay with an 
area law. 
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Figure 3. A comparison between square Wilson 
loops and U-shaped operators of the same size at 
(ftj, fa) = (0.75245, 0.16683) where TZ b yfa ~ 27. 



Measuring (W(r, t)) and (C r ) yields an esti- 
mate of the string breaking scale TZb- While (W), 
as it turns out, does not depend on TZb at the 
probed scales, the U-shaped operator is very sen- 
sitive to it: (U(r,t)) drops off rapidly as TZb in- 
creases (see Fig|3J|. This is not unexpected: when 
TZb — > 00 the theory becomes a pure gauge theory, 
therefore (U) = 0. 

In order to investigate in more detail such a 
dependence on the string breaking scale it is con- 
venient to study a new family of gauge-invariant 
operators, that we call almost closed Wilson lines 
G r (l,t), which somehow interpolate between the 
U-shaped ones and the Wilson loops (see Fig0}. 
It turns out that the vacuum expectation value 
(vev) of these operators does not depend very 
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Figure 4. The almost closed Wilson lines interpo- 
late between the Wilson loops and the U-shaped 
operators. 



much on the position of the pair of matter sources 
• • along 7. This can be expressed with the 
graphical equation 



(7) 



(W(l,l))/(G r (l,l)) 
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depend, within the errors, on the size I, as FigEl 
demonstrates. Although this numerical fact can- 
not directly promoted to the rank of a physical 
property of the continuum theory because of un- 
balanced self-energy divergences between numer- 
ator and denominator, by combining the almost 
closed Wilson line G r (l,l) with the straight line 
operator C r one is led to consider the following 
gauge- invariant, UV finite quantity 



(C r ) (G r (lJ)) 
(W(l,l)) 



f(r) 



(8) 



which, according to the results of FigEI seems 
to depend only on r. The behaviour of such a 
function is drawn in Fig|5] More insight into f(r) 
will be gained in the next section. 
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' Figure 6. The function /(r) defined in Eq.lJSJ 

evaluated at (/3 G ,Pi) = (0.735,0.224). The con- 
tinuous line is a fit to j,eT mT ' . The parameters c 

Figure 5. Ratio between a square Wilson loop and m wil1 be relatcd to othcr Physical quantities 

(W(l, I)) and an almost closed square Wilson line °^ tne system. 

of the same size (G r (l,l)). One clearly sees that 

such a ratio does not depend on the size I, but 

only on the distance r between the end points. 

3. FACTORIZATION 

Evaluating the ratio {W(l, l))/(G r (l, I)} one There is no reason to believe that the almost 

discovers the surprising result that it does not closed Wilson lines play some special role in gauge 
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Figure 7. Splitting a Wilson loop into two open 
Wilson lines. 



theories, thus Eq.© should be regarded as re- 
flecting a more general property of whatever open 
Wilson line. An obvious generalisation which 
comes to mind is to consider an arbitrary planar 
Wilson loop W(7) (see FigEJ, cut the closed path 
7 in two points placed at a distance r and satu- 
rate the end points of the resulting open paths a 
and (3 with matter fields. In this way one is led 
to write down the following UV finite quantity, 
generalising Eq. (JSJl 



(GV(a)> (G r {f3)) 



■0)) 



f(r) 



(9) 



In order to lend some support to such a proposal 
we measured the ratios (C r ) (U(r,r)) / (W(r,r)) , 
(U(r,r))y(W(2r,r)) and (G^l) ) 2 /{W(r, r)>, 
where L denotes the path obtained by cutting a 
square in half along its diagonal. Some results arc 
presented in Fig|HJ 

Though these numerical tests of Eq.© arc 
manifestly very good, there are regions of the pa- 
rameter space where the agreement is less good. 
The general rule seems to be that Eq.© is ful- 
filled, at least approximately, whenever the size 
of the involved loop is less than the string break- 
ing scale. But even in this region Eq.© cannot 
reflect an exact property of the gauge theories. 
The reason is very simple: the vev of the Wilson 
loop W{"f) in the denominator depends explicitly 
on the shape of 7 and this dependence cannot be 
counterbalanced by the shape effects of the nu- 
merator, being the splitting 7 — > a+f3 completely 
arbitrary. Such a fact will be strengthened in 
the next section by an explicit calculation of the 
universal shape effects generated by the quantum 
fluctuations of the underlying confining string in 
the U-shaped open Wilson line. 
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Figure 8. Testing Eq.©. The CPs correspond 
to {G^2 r (Z.)) 2 , the *'s correspond to (U(r,r)) 2 
and the x's to (C r ) (U(r,r)), finally the +'s and 
the continuous line have the same meaning as in 
FigEl but here the data are taken at (Pg^Pi) = 
(0.75245,0.16683). Data corresponding to the 
same value of r are displaced for clarity. 



In conclusion, Eq.© has no chance to be exact. 
Rather, it should be viewed as the first term of 
an expansion in some scale parameter, perhaps 
1/aAy. Note that if 7 is scaled up to 00 keeping 
its shape fixed one gets (a generalisation of) the 
Frcdcnhagcn-Marcu order parameter [7] 



lim 

r — >oo 



(G r (a)) (Gr(0)) 
(W( 7 )> 



fixed 
shape 



(10) 



which vanishes in the Coulomb-like phase while 
it is different from zero in the confinement /Higgs 
phase. 

Eq. © suggests the first correction to p should 
be a shape-independent term. Being only a func- 
tion of the distance r of the two end points where 
are sitting the matter fields, it is conceivable that 
f(r) has the asymptotic expansion of a two-point 
correlation function, namely 



„d-2 



+ 



(11) 
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which fits well to the numerical data as FigEJand 
[H] show. The parameter p is too small and can 
be put safely equal to zero in these fits. A better 
evaluation of p will be obtained shortly. 

It has to be noted that the "mass" m, though 
is a physical quantity needs not to correspond to 
the mass of any state of the physical spectrum, 
being associated not to a single correlator, but 
to a suitable rational function of three of them. 
A previously unsuspected relation between m and 
the string breaking scale will be uncovered below. 

Neglecting for simplicity the pre-exponential 
factor l/x d ~ 2 and using as input the Ansatz (0 
we get 

(G r (7)) - (Wm(c e ^" m ) r + p)/d . (12) 

The first term shows that open Wilson lines obey 
a perimeter-area law which differs from that of 
the ordinary Wilson loops in two respects. First, 
the proportionality constant is different. Sec- 
ondly, there is no analogue of the contribution 
of that part of the boundary where the string is 
free, encoded in the coefficient Cf of r in the ex- 
ponential, which clearly is c/ = p — A — m. 

In order to find the link between m and IZi, 
we alluded to above, the clue is provided by the 
observation that Eq.JSJl implies the consistency 
condition 

(U{r,t)) 2 ^ (U(r,2t)) (£,) (13) 
or, diagrammatically, 

(qq^d ■)<]) . (w) 

It gives at once 

cf m , . 

c = — ; Cf = — — ; m = crTL b . (15) 

Cu ^ 

Combining Eq.s l|12|) an( l G5t we get finally 
Eq.(@J, which was the major purpose of this talk. 

In a nutshell, we may conclude that the sur- 
prising factorization properties of the vev of the 
open Wilson lines tell us nothing more than their 
behaviour at intermediate scales is dominated by 
the exponential law (J2J. 

One could try to look for a more detailed so- 
lution of (JHJ by taking into account the pre- 
exponential factor \/ r d ~ 2 of Eq. (|llfl . Note how- 
ever that such a kind of correction is of the 
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Figure 9. The U-shaped Wilson line at (/3 G , (ii) = 
(0.735,0.224). The change in slope is due to the 
non- vanishing of p, which is the only free param- 
eter of the fit. 



same order of the already mentioned shape ef- 
fects, which we know do not obey factorization, 
thus using factorization to find the sub-dominant 
corrections to the exponential behaviour is no 
longer justified. 

An interesting property of the first two terms 
of Eq. is that all the parameters but one (the 
Fredcnhagcn- Marcu order parameter) can be de- 
termined by evaluating the vev of ordinary Wil- 
son loops or straight Wilson lines, hence the fits 
to this equations are one-parameter fits. As an 
example a fit to the U-shaped Wilson line is pre- 
sented in FigEl Note that the dependence on p is 
strongly enhanced by an exponentially increasing 
factor. 

4. UNIVERSAL SHAPE EFFECTS 

It is known that the description of the vev of 
Wilson loops as a simple exponential like in |JTJ 
is not sufficiently accurate. Even in gauge theo- 
ries coupled with matter the confining string, as 
long as its length does not exceed TZb, can vi- 
brate freely exactly like in the absence of dynam- 
ical matter |21 1 j . It is then natural to conjee- 
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ture that also in the open Wilson line G r (a) one 
could detect universal shape effects generated by 
the quantum fluctuations of the confining string 
like in the ordinary Wilson loop. The important 
difference is that now the boundary conditions 
of the string are fixed only along the open path 
a, whereas are free in the part of the boundary 
formed by the straight line connecting its end- 
points. 

In the case of the U-shaped operator it is easy 
to evaluate the contribution of these fluctuations. 
To be definite, consider the operator U(r, £), cor- 
responding to the case r = I of Fig0J The bound- 
ary conditions are chosen accordingly. The uni- 
versal shape effects are encoded in the quantity 

una 



[det(-c> 2 )] = \-d 2 \-^ , (16) 

where d 2 is the Laplacian of the 2d bosonic field 
describing the string displacements in the d — 2 
transverse directions with respect to the minimal 
surface. We have 

i- 92 i=nn^ 2 +^0 2 ]> m 

m>0 n>0 

which clearly needs rcgularization. A simple 
method is based on the assumption that reg- 
ularised scries (or products) fulfil some formal 
properties of the absolutely convergent series. In 
particular, we assume the two relationships 



n>0 



2*[J> 

n>0 



n>0 



n>0 



a-1 
n=l 



n>0 



(18) 



(19) 



where the apex indicates the regularised expres- 
sion and the finite sum X)n=i = Pk+i ( a ) i s a 
polynomial of degree k + 1. For our purpose we 
need only to know Pi (a) = a — 1 and Pzfa) = 
a(a — l)/2. Assuming they are valid for any real 
a we get at once 



^n fe = 2 fe P fe+1 (i)/(2 fc+1 -l), 



(20) 



n>0 



which turns out to coincide with the value of the 
Riemann zeta function £(— k). In this spirit we 



get the identities 

' ' 1 

[]c = exp[logc ]Tn°] = — (21) 

ti>0 n>0 * 

' ' 1 

H c = exp[logc ^(n--)°] = l (22) 

' ' 1 

JJn= V2^ , l[(n--) = V2 (23) 

n>0 n>0 

Which can directly applied to Eq. l|17|l . yielding 



1 + ¥ 



i-^nn 

m>0n>0 . 

Using the known identity 



na „ 
+ ^ 



n>0 



lira 



2na 



(24) 



.^_ e «hraj ^25) 



combined with Eq.s (0,1(211) and leads us 
finally to the finite result 



^ = ^9-* IK 1 



q 



), q = e 



(26) 



m>0 



which can be written in a number of equivalent 
ways: 



i) 2 



h(r) _ rj(r/2) _ r,(-2/r) 



2t?(t) y/2 v (T) »?(-l/r) 



(27) 



where r = i~, 9± is a Jacobi ^-function and r/(r) = 

riri>o(l — q n ) 1S the Dedekind's eta function. 
In conclusion, the first term of the asymptotic 
expansion |@J applied to the U-shaped Wilson line 
can be written in the form 

(U(r, t)) oc F(t/r) e - rt "-( 2 *+<-) >-rn b «/2 + . . . ^ 

As a check, it is easy to see that the static po- 
tential defined as V(r) = — lim^oo log((U(r, t))) 
contains the 1/r Liischer term with the right co- 
efficient. 

Observing these effects is very challenging. 
Strange though it may seem, boundary terms con- 
stitute the major obstruction to detect universal 
string signals. 
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Indeed finite size effects in interfaces, where the 
involved string world-sheet has no boundary, is 
the first place where these effects were clearly seen 

G3 

Similarly, in the Polyakov loop correlators the 
boundaries contribution does not depend on their 
distance; indeed in such a system are concen- 
trated the main theoretical efforts to understand 
the physics of the confining string, after the in- 
troduction of the variance reduction method |18| . 
For a recent discussion on this argument and a 
complete list of references see [T§] . 

In the case of Wilson loops the boundary term 
increases with the loop size, hence it is much more 
difficult to achieve accurate results. As a matter 
of fact only in very simple systems, like in 3D Z2 
gauge model |20| or even in percolation processes 
|21| high precision estimates have been reached. 

In the case of open Wilson lines the situation 
dramatically worsen, the reason being that there 
are two competing requirements: the shape ef- 
fects would be visible in sufficiently large open 
Wilson lines provided the string does not break, 
hence one has to chose a region where IZb is large, 
however this implies, as Ea. (|28|) shows, an expo- 
nential weakening of the signal, because of the 
contribution of the free boundary. Studies are 
under way to gain more insight into these opera- 
tors. 
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